Dispersive force between dissimilar materials: geometrical effects 
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We calculate the Casimir force or dispersive van der Waals force between a spherical nanoparticle 
and a planar substrate, both with arbitrary dielectric properties. We show that the force between a 
sphere and a plane can be calculated through the interacting surface plasmons of the bodies. Using 
a Spectral Representation formalism, we show that the force of a sphere made of a material A and a 
plane made of a material B, differ from the case when the sphere is made of B, and the plane is made 
of A. We found that the difference depends on the plasma frequency of the materials, the geometry, 
and the distance of separation between sphere and plane. The differences show the importance of 
the geometry, and make evident the necessity of realistic descriptions of the sphere-plane system 
beyond the Derjaguin Approximation or Proximity Theorem Approximation. 



I. INTRODUCTION 

The origin of dispersive forces like Casimini and van 
der Waals (vW) forces^* 3 -, between atoms and macro- 
scopic bodies may be attributed to electromagnetic in- 
teractions between their charge distributions induced by 
quantum vacuum fluctuations, even when they are elec- 
trically neutraiiiSi&i In a first approximation, the charge 
distribution of neutral particles may be represented by 
electric dipoles. This dipole approximation was em- 
ployed by London to calculate the non-retarded van der 
Waals interaction potential V v yj(z) between two identi- 
cal polarizable molecules by using perturbation theory 
in quantum mechanics 5 . Later, Casimir studied a sim- 
pler problem^: the force between two parallel conduct- 
ing plates separated by a distance z, due to the change 
of the zero-point energy of the classical electromagnetic 
modes. He found an interaction energy per unit area, 
y A (z) = - (n 2 he / 720) (1/z 3 ), where c is the speed of light. 
In 1956, Lifshita 3 . extended the theory of Casimir to di- 
electric semi-infinite slabs. A decade later, it was shown 
that the Lifshitz's formula could be rederived from the 
zero-point energy of the surface modes of the slaba£»Z»£. 

In recent years, accurate experiments&iSiiiii^ii 3 . have 
been performed to measure the Casimir force which is 
one of the macroscopic manifestations of the fluctua- 
tions of the quantum vacuumi. Most of these experi- 
nients 10 ' 11,12 ' 13 have been done using a spherical surface 
and a plane, instead of two parallel planes, as originally 
proposed Casimir 1 himself. The interpretation of the 
Casimir force in these experiments is based on the Prox- 
imity Theorem which was developed by Derjaguin and 
collaborators^ to estimate the force between two curved 
surfaces of radii R\ and i?2- The Proximity Theorem or 
Derjaguin Approximation (DA) assumes that the force 
on a small area of one curved surface is due to locally 
"flat" portions on the other curved surface, such that, 
the force per unit area is 



T{z) = 2tt 



R1R2 

Rl + i?2 



V(z), 



when Ri = R, and R2 — > 00, the problem reduces to the 
case of a sphere of radius R and a flat plane, that yields 
to T{z) = 2ttRV(z). The force obtained using the DA 
is a power law function of z, where at "large" distances 
is proportional to z~ 3 , while at short distances T(z) oc 
z~ 2 . This theorem is supposed to hold when z <C R\, R2, 
however, the validity of the DA has not been proved yet. 

In this direction, Johansonn and ApelU 5 . studied the 
vW force between a sphere and a semi-infinite slab. They 
calculated the electromagnetic stress tensor associated 
to the electric field correlation of the system from the 
Green's function of the problem expressed in bispherical 
coordinates. They concluded that for small separations 
the behavior of the attractive force is consistent with the 
DA. However, one drawback of this formalism is that in 
bispherical coordinates the section surfaces become pla- 
nar (a point) for small (large) values of 2, so that, ar- 
bitrary values of ratio of z/R cannot be considered to 
study the sphere-plane system. 

Casimir showed that the energy li(z), between two par- 
allel perfect conductor planes can be found through the 
change of the zero-point energy of the classical electro- 
magnetic field—, like 



U ^ = 9 X^ 2<z ) ~ Ui ^ z ~* °°)]' 



(1) 



where u>i(z) are the proper modes that satisfy the bound- 
ary conditions of the electromagnetic field at the planes 
which are separated a distance z. For real materials, Lif- 
shitz obtained a formula to calculate the Casimir force 
between parallel dielectric planes^. The force per unit 
area, according to the Lifshitz formula, is 
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where V(z) is the Casimir energy per unit area between 
parallel planes separated by a distance z. In the limit, 



with £ Q = [rfr% exp(2ifcz)] _1 , where r" is the reflection 
amplitude coefficient of the plane j for the electromag- 
netic field with a-polarization. Here, j = 1, or 2, a — s 
or p, q — uj/c, u> is the frequency, c is the speed of light, 
q = (Q, k) is the vacuum wavevector with projections 
parallel to the surface Q, and normal to the surface k, 



k — k + iO + , and g is the photon occupation number of 
the state k. 

The Lifshitz formula depends only on the reflection co- 
efficients of the planes and the separation between them. 
This reflection coefficients can be calculated for arbitrary 
material using the surface impedance formalism for semi- 
infinite slabsi^iLi^ and for finite slabs using, for exam- 
ple, the transfer matrix formalismi2i2£. In addition, one 
observes that the Lifshitz formula is symmetric under 
the change of index j = 1 — ► 2, it means that the force 
is indistinguishable under the change of planes which is 
natural given the symmetry of the system. 

The force between a sphere and a plane made of ar- 
bitrary dielectric materials is commonly calculated using 
the DA, where the Casimir energy V{z) is found using the 
Lifshitz formula. As a consequence, the force between a 
sphere and a plane is indistinguishable if the sphere is 
made of a material A and the plane is made of B, or 
if the sphere is made of B and the plane is made of A. 
Furthermore, the DA assumes that the interaction of a 
sphere and a plane can be calculated as the interaction 
of a plane, and a set of planes, and evaluates the force 
by adding contributions of various distances, as if the in- 
teraction between planes were independent. However, it 
is well known that the Casimir force is not an additive 
quantityiSi. 

In this paper, we show that the geometry of the system 
is important. We employ a method based on the deter- 
mination of the proper frequencies of the system based 
on a spectral representation formalism^. By calculating 
the zero-point energy of all the interacting surface plas- 
mons of a sphere-plane configuration, we obtain that the 
force between dissimilar materials differs if the sphere is 
made of a material A and the plane is made of B, and 
vice versa. This provides an insight into the range of va- 
lidity of the DA or Proximity Theorem. Also, we find 
a general formula to calculate the Casimir force differ- 
ences between arbitrary materials when the bodies are 
separated at least one radii of the sphere. In particular, 
we study the specific case of aluminum and silver. We 
restrict ourself to the nonretarded limit, this means, to 
the case of small particles at small distances. 



In the case of two parallel planes, in the absence of 
retardation, the proper electromagnetic modes satisfy the 
following expression 6 



e(w) + 1 



1 2 



e(«) - I 



exp[2A:z] -1 = 



(2) 



for any value of sj k sj oo. The roots u>i(k) are identified 
as the oscillation frequencies of the surface plasmons. To 
find the proper modes, it is necessary to choose a model 
for the dielectric function of the planes. To illustrate the 
procedure for metallic planes, let us employ the plasma 
model for the dielectric function, as 



eM = !--§, 



(3) 



with Up the plasma frequency. Substituting Eq. J3J in 
Eq. (J2J, we find two proper modes that, at large distances 
(kd » 1), are given by 



w± « ^ (l ± \ exp[-fcd] - ~ exp[-2fcd] ± 



(4) 



and in the limit kd — ► oo, we obtain u± = ui v j \/2 which 
are the frequencies of the surface plasmon of each half 
space. The zero-point energy per unit area will be given 
by 



£o(z) 



5^jfl>(*)**' 



(5) 



Now, taking the difference of the zero-point energy when 
the planes are at a distance z, and when they are at 
infinite, we obtain the interaction energy of the system 
per unit area, like 



V(z) 



hw„ 



y/2 IGz 2 ' 



(6) 



which depends on the energy of the surface plasmons of 
the half spaces and the distance between planes. 



II. ZERO-POINT ENERGY OF THE 
INTERACTING SURFACE PLASMONS 

A. Parallel planes 

In 1968, van Kampen and collaborators 6 showed that 
the Lifshitz formula in the nonretarded limit is obtained 
from the zero-point energy of the interacting surface plas- 
mons of the planes. Later, Gerlachi did an extension 
showing that also in the retarded limit the Lifshitz for- 
mula is obtained from the zero-point energy of the inter- 
acting surface plasmons. A comprehensive derivation of 
the Lifshitz formula using the surface modes is also found 
in Ref. 



B. Sphere and plane: dipolar approximation 

Now, let us consider the case of a sphere of radius a 
and dielectric function e s (w), such that, the sphere is 
located at a minimum distance z from the plane of di- 
electric function e p (u>), as shown in Fig. 1. The quantum 
vacuum fluctuations will induce a charge distribution on 
the sphere which also induces a charge distribution in the 
plane. Then, the induced ^m-th multipolar moment on 
the sphere is given by22 



hmW = <*, m (w) [tc» + ^»] , 



(7) 



where V^ c (uj) is the field associated to the quantum vac- 
uum fluctuations at the zero-point energy, VJ™ b (w) is the 




FIG. 1: Schematic model of the sphere-plane system. 



induced field due to the presence of the plane, and ai m (iv) 
is the Zm-th polarizability of the sphere. 

To illustrate the procedure in detail, first we work in 
the dipolar approximation, i.e., when 1 = 1. In section 
III, we will calculate the proper modes taking into ac- 
count all the high-multipolar charge distributions. Using 
the method of images, the relation between the dipole 
moment on the the sphere p s (u>), and the induced dipole 
moment on the plane p p (uj), is 



P P M 



1 - £ p(^) 1 
1 + e p {u>) 



■Ps(w)- 



(8) 



Here, M = (—1,-1,1) is a diagonal matrix whose ele- 
ments depend on the choice of the coordinate system. 
From Eq. (JZJ), one finds that the total induced dipole 
moment on the sphere is 



fi(w) = o(w) £ va » + T-;p P H 



(9) 



where the interaction between dipoles is coupled by the 
dipole-dipole interaction tensor 

T= (3rr-r 2 I)/r 5 . 

Here, I is the identity matrix, and r is the vector be- 
tween the centers of the sphere and the dipole-charge 
distribution on the plane. From Fig. 1, we find that r — 
(0,0,2(z + a)), such that, M-T= (-1/r 3 , -l/r 3 ,-2/r 3 ) 
is a diagonal matrix. Substituting Eq.© in Eq.@, one 
finds 



a(u>) 



I + / C (u/)M • T 



(10) 



where G(u>) is a diagonal matrix and 

/c(«) = [l-e p («)]/[l + e p («)]. 

Multiplying Eq. (|10|l by a 3 , one finds that each term of 
&(oj) is dimensionless and has two parts: the left-hand 
side which is only associated to the material properties of 



the sphere through its polarizability l/a(u>) = a 3 /a(to), 
and the right-hand side which is related to the geometri- 
cal properties of the system trough a, and z. The right- 
hand side also depends on the material properties of the 
plane trough the function / c (w). 

The eigenfrequencies of the sphere-substrate system 
must satisfy Eq. (|10J) . and are independent of the ex- 
citing field, in this case E vac (cu). These frequencies can 
be obtained when the determinant of the matrix in the 
left-side of Eq. (|10|l is equal to zero, i.e., detG(w) = 0. 
In this particular case, the eigenfrequencies can be found 
from the following expression 



1 



/c(w)o 3 



l 2 



a(w) [2{z + a)f 



1 



2f c (u>)a 3 



= 0. (11) 



_a(u) [2(z + a] 

At this point, it is necessary to consider a model for the 
dielectric function of the sphere to find the proper modes. 
Again, we illustrate the procedure for metallic sphere and 
plane, using the plasma model of Eq. ©. Moreover, we 
use the fact that the polarizability of the sphere, within 
the dipolar approximation and in the quasi-static limit, 
is given by 



a{io) 



M - 1 



£ s (w) + 2 

In Fig.Et, we plot, l/a(w) and / c (w)[2(l + z/a)]~ 3 as 
a function of lu/lo p for different values of z/a. The proper 
frequencies of the system are given when the black-solid 
line and the color-dotted lines intersect each other. We 
observe that there are two different proper modes for each 
term on Eq. I|ll(l . giving a total of six modes for each 
z/a. The proper modes at the left-hand side in Fig. [SJi 
correspond to the surface plasmon of the sphere that we 
denote with cj+. These modes are red-shift as the sphere 
approaches the plane because of the interaction between 
them. When the separation z/a increases, the modes w+ 
go to cjp/V3, which is the surface plasmon frequency of 
an isolated sphere. On the other hand, the proper modes 
at the right-hand side correspond to the surface plasmon 
of the plane and we denote them with a;_ . These modes 
are blue-shift as the sphere approaches the plane, also 
because of the interaction between the sphere and the 
plane. In this case, when the separation z/a increases, 
the mode w_ goes to uj p /y2 which is the frequency of 
thesurface plasmon of the isolated plane. This behavior 
of the modes is clearly observed in Fig.[2), where lu + /uj p 
and lo-/lu p are plotted as a function of z/a. This shows 
that the the zero-point energy of the sphere-plane system 
is directly associated to the interacting surface plasmons 
of the sphere and the plane. 



III. EXACT CALCULATION OF THE 

ZERO-POINT ENERGY OF A SPHERE ABOVE 

A PLANE 

In this section we calculate the proper electromag- 
netic modes of the sphere-plane system, within the non- 
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FIG. 2: (Color online) (a) Plot of l/5(w) in solid line and 
/ c (u>)[2(l + z/a)]~ :i as a function of w/w p for different values 
of z/a. (b) Proper modes as a function of z/a. 



retarded limit, by including all the high- multipolar in- 
teractions. These proper modes are calculated using a 
Spectral Representation formalism22i25i2& ) then, we cal- 
culate the zero-point energy using Eq. QJ. The Spectral 
Representation (SR) formalism uses the fact that the con- 
tributions of the dielectric properties of the sphere can be 
separated from the contributions of its geometrical prop- 
erties. The latter allows to perform a systematic study 
of the sphere-plane system that may be very helpful for 
experimentalists. The details of the SR formalism can be 
found in Refs. 122 andl27l here we only explain it briefly. 
We now calculate the proper modes of the system in- 
cluding all the high-multipolar charge distributions due 
to quantum vacuum fluctuations. Using the method of 
images, we find from Eq. that the Im-th multipolar 
moment on the sphere is 



E 

I'm' 



4ir6a'$n 



l m 



(2J + 1)ojw(w; 



+ f c (u)A\™ (z) 



tETM, (12) 



where A\™ (z) is the matrix that couples the interac- 
tion between sphere and substrate^. The expression for 



the interaction matrix A\™ (z) is given in the Appendix 
section, and it shows that A\™ (z) depends only on the 
geometrical properties of the system. If we have a ho- 
mogeneous sphere, its polarizabilities are independent of 
the index to, and are given by2£ 



ai(w) 



l[e s (tu) - 1] 
1[e s (uj) + 1} + 1 



,21+1 



nio 



nio - u{u) 



,21+1 



(13) 



where ni — 1/(21 + 1), and u(lo) = [1 — e 8 (w)] -1 . No- 
tice that in the above equation for the polarizability, the 
dielectric properties of the sphere are separated from its 
geometrical properties. By substituting the right-hand 
side of the above equation, we rewrite Eq. (|12|l as: 

E { - «m^<w + <rV ^ ^fy/2 = 



I'm, 



21+1U/2 



(la 2i+1 ) 

47T 



-^ C M, (14) 



where the lin, I'm'-th element of the matrix H(w, z) is 

H l lr ^ (w,z) = ni>o5wS mm > + 

f c (u;)^ip.(ll')^A l t'(z). (15) 

The proper modes are given when the determinant of the 
quantity between parenthesis in Eq. (|14fl is equal to zero, 
det [-u(u>)I + H(w, z)] = 0, that is 



G{u, z) = JJ[-u(w) + n,(w, z)] = 0, 



(16) 



where n s (u>, z) are the eigenvalues of M(u>, z). 

The Spectral Representation formalism can be applied 
in a very simple way. First, we have to choose a substrate 
to calculate the contrast factor / c (w), with this, we con- 
struct the matrix H(w, z), and we diagonalize it to find 
its eigenvalues n s (cj, z). We can repeat these steps for a 
set of different values of z. Notice that the eigenvalues 
of H(w, z) can be found without to do any assumption 
on the dielectric function of the sphere. Once we have 
the eigenvalues as a function of z, we have to consider 
an explicit form of the dielectric function of the sphere, 
such that, we calculate the proper electromagnetic modes 
U) s (z) trough the relation u(lo s ) — n s (uj,z). 

To illustrate the procedure, we use again the plasma 
model for the dielectric function of the sphere, and a 
plane with a constant dielectric function. In this case, 
the proper modes are given by 



u g (z) = u) p \/n 3 Jz), 



(17) 



with s = (l,m), and where the eigenvalues n 3 (z), are 
independent of the frequency. The zero-point energy, ac- 
cording with Eq. Q , is given by 



fko 



£( z ) = -^E [\AW* 



'n m 



(18) 



l.ni 



where hujp^/nio 



hjjJ v \Jlj(2l + 1), are the energies asso- 
ciated to the multipolar surface plasmons resonances of 
the isolated sphere. Here, I = 1, 2, . . . , L, where L is the 
largest order in the multipolar expansion. For I = 1 we 
have that huj p ^/nio = fiw p /vo, that is the surface plas- 
mon of the sphere in the dipolar approximation. When 
I — > oo, we have that huj p ^nio ~* fiw p /v2- 

Now, in the presence of the dielectric plane, the proper 
modes of the isolated metallic-sphere are modified ac- 
cording to Eq. I|15l) . These proper modes are red-shift 
always as the sphere approaches the plane, and this shift 
depends on the separation z/a. We find that the inter- 
action energy is proportional to [1/2(1 + z/a)]' 3 , where 
j3 = I + I' + 1 is an integer between 3 and 2L + 1. As 
z/a — > 0, more and more multipolar interactions must 
be taken into account and f3 — > 2L + 1. Finally, when 
the sphere is touching the plane, one would expect that 
L — > oo, so that also (3 — > oo, and the energy at z = is 
proportional to (1/2) 7 , with 7 — > 00, it means the energy 
becomes also infinite. 

It is possible to find the behavior of the energy as a 
function of z/a for any plasma's sphere independently of 
its plasma frequency, like £ = £/hu> p . Furthermore, the 
force can be also studied independently of the radius of 
the sphere, since £{z/ a), one finds that 



F = - 



d£{z/a) d£(z/a) d{z/a) 



Oz 



d(z/a) dz 



(19) 



and then, we can define also a dimensionless force like 
aF = aF/huip. In summary, we have shown that given 
the dielectric properties of the substrate the behavior of 
the energy and force is quite general, showing the poten- 
tiality of the Spectral Representation formalism^. 

In Fig. [3J we show the dimensionless energy, £, and 
force aF, for a sphere and a plane both described by 
the plasma model with the same plasma frequency. We 
show £ and aF as a function of z/a when all multipolar 
interactions are taken into account, as well as only dipo- 
lar, and up to quadrupolar interactions are considered. 
The largest multipolar moment to achieve convergence 
of the energy at a minimum separation of z/a = 0.1 
is L = 80. We observe that £ and aF are power law 
functions of (z/a) - ' 3 , where j3 depends on the order of 
the multipolar interaction considered. Only within the 
dipolar approximation, £ and Fa behaves like (z/a)~ 3 , 
and (z/a) -4 for any z, respectively. Let us first ana- 
lyze the system when only up to quadrupolar interac- 
tions are taken into account. In that case, the energy 
as well as the force, show three different regions: (i) at 
large distances, z > 5a, only dipolar interactions are im- 
portant, (ii) when 5a > z > 2a, we found that dipolar- 
quadrupolar interactions become important and £ and 
Fa behaves like (z/a)^ 4 , and (z/a)~ 5 ; while (iii) at small 
distances z < 2a the quadrupolar-quadrupolar interac- 
tion dominates, and £ and Fa behaves like (z/a)~ 5 , and 
(z/a)~ 6 , respectively 

However, as the sphere approaches the substrate the 
interaction between high-multipolar moments becomes 
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FIG. 3: (Color online) (a) Energy £, and (b) force aF as 
a function of z/a, with L = 80 (dotted line), L = 2 (solid 
line), and L = 1 (dashed line), (c) Force difference, \(F LH — 
F LW )/F LH \, using multipolar moments LH = 80 and LW = 
2 (solid line), LH = 80 and LW = 1 (dashed line), finally 
LH = 2 and LW = 1 (dot-dot-dashed line). 



more and more important. The attractive force suddenly 
increases more than four orders of magnitude, as compare 
with the dipolar and quadrupolar approximations. In 
Fig. I3J; we show the difference of the dimensionless force, 
\(pLH _ F LW )/F LH \, as a function of the separation, 
with LH and LW the highest and lowest-multipolar mo- 
ments taken into account. We observe that at distances 
z > 2a the force can be obtained if up to quadrupo- 
lar interactions are considered. We also found that for 
z > la the interaction between the sphere and the sub- 
strate can be modeled using the dipolar approximation 
only, like in the Casimir and Polder models. However, 
at separations smaller than z < 2a, the quadrupolar ap- 
proximation also fails, and it is necessary to include high- 
multipolar contributions. For example, at z = a the 
quadrupolar approximation gives an error of about 10%, 



while the dipolar approximation gives an error of about 
40%. At smaller distances like z = a/2, the quadrupo- 
lar approximation gives an error of about 40%, while in 
the dipolar approximation the error is 70%. At z = 0.1a 
both approximations give an error larger than 90%. To 
directly compare the dipolar and quadrupolar approxi- 
mations, we also plot in Fig. [3: the force difference when 
LH = 2 and LW = 1. In this curve, we clearly observe 
that even at z = 5a there are differences between the 
dipolar and quadrupolar approximations of about 5%. 

The large increment of the force at small separations 
due to high-multipolar effects could explain the physi- 
cal origin of, for example, the large deviations observed 
in the deflection of atomic beams by metallic surfaces 29 , 
as well as some instabilities detected in micro and nano 
devices. Furthermore, to compare experimental data of 



some experiment: 



,12,30 



with the DA, it has been needed to 



make a significant modification to the Casimir force. In 
Ref. [13 the authors found a lar ger f orce than the one cal- 
culated using DA, while in Ref. yfjthe authors measured 
a smaller one. Both groups argue that the discrepan- 
cies between theory and experiment are due to the same 
effect: the roughness of the surface. However, the dis- 
cussion around this point has not sense until the DA 
be proven. Indeed, this ambiguity on the interpretation 
of the measurements indicates that the force between a 
spherical particle and a planar substrate involves more 
complicated interactions than the simple dipole model 2 ^ 
or the DAii. For example, such deviations might be 
attributed also to high-multipolar effects, however, new 
experiments have to be perform to prove the later. 

In this paper, we propose a way to elucidate experi- 
mentally two important issues to the Casimir force: (i) 
the relevance of the geometry via the surface plasmon 
interactions between bodies, and (ii) the relevance of the 
high-multipolar interactions in the sphere-plane model. 
In the next section, we study the case of the casimir force 
between a sphere and a plane of dissimilar materials that 
can help us to clarify these important issues. 



IV. CASIMIR FORCE BETWEEN DISSIMILAR 
MATERIALS 

Using the formalism presented in the previous section, 
we calculate the Casimir force for a sphere and a plane 
made of dissimilar materials. To illustrate our results, 
we choose silver and aluminum whose dielectric function 
is described by the Drude model 



e(w) 



lu(lu + i/r) ' 



(20) 



with plasma frequencies hui p = 9.6 eV, and 15.8 eV, 
and relaxation times (rtUp)" 1 = 0.0019, and 0.04, respec- 
tively. In this case, we employ the Argument Principle 
(AP) to calculate the force. The AP was first used by van 
Kampen 6 and Gerlach 7 to calculate the force between 



planes. The reader can consult Ref. 4 for a comprehen- 
sive explanation of the use of the AP within the context 
of Casimir forces. 

Then, employing the Argument Principle method we 
find that the sum over proper modes is given as 



*(*)=E 



hw s (z) h f 8G(lu,z) 1 
2 = 4^/ c W duj G(u>,z) 



-' : rdiU , 

(21) 



where G(cu, z) is defined in Eq. I|16|) . such that 
£(z) = — ^T j> w— log [-• u(w) + n,(w, z)]dw. (22) 

s " C 

The closed contour C is defined as the part along the 
imaginary axis plus the part along the semicircle on the 
positive real axis. In the limit for a infinite semicircle 
radius, the integral along the semicircle does not depend 
on the separation z, therefore, it does not contribute to 
the force. Then, it follows that the interaction energy of 
Eq. 0J is given by 



U{z) 



Am 



E 



dw \og[— u(uj) + n s (oj,z)], (23) 



where in the last step we have performed a partial inte- 
gration. The force is found as in Eq. I|19|l . where we first 
perform the derivate with respect to z of Eq. (|23|l and 
then the integral with respect to ui. 

In Fig. 01 we show the calculated force for an aluminum 
sphere of arbitrary radius above a silver plane, ai*Al/Ag> 
in solid line, and the same for a gold sphere above an 
aluminum plane, ai^Ag/Ab i n dashed line. We find that 
the difference is small and cannot be directly observed. 
Then, in Fig. 2t> we show the difference of the force nor- 
malized with the average as a function of z/a, given by 



AJ" = 2 



-Pa/b - Fb/a 



■"A/B 



i 7 ! 



B/A 



(24) 



with A=A1, and B=Ag. Here, we observe that the dif- 
ference of force varies as a function of z/a, such that at 
small distances, (z < 0.5a) the difference under the inter- 
changed of materials between the sphere and the plane, 
is less than 3%. Then, the difference increases with the 
minimum separation, and at z > a, it has almost reached 
their maximum. For these particular materials, the max- 
imum is about 4.86%, independently of the radius of the 
sphere. We find that this value of A3 7 is a consequence 
of the geometry of the system, as we explain below. 

In the previous section, we show that for a minimum 
separation larger than 2a, the force can be described with 
quadrupolar interactions between the sphere and plane. 
This means that for z > 2a the force varies with the 
same power-law function. This conclusion can be easily 
deduce from Fig. |2t. We also show in section II, that 
the Casimir force is given by the interaction of the sur- 
face plasmons of the sphere and plane. In Fig. [Sb, we 
observe that the proper modes of the system for z > a 




FIG. 4: (Color online) (a) Force as a function of z/a, for 
a sphere made of Al over an Ag plane (solid line), and for 
a sphere made of Ag over an Al plane (dashed line), (b) 
Difference of the force between a sphere made of Al over an 
Ag plane, and the force for a sphere made of Ag over an Al 
plane. 



are approximately given by hw p /^/2 and hu)p/V3, when 
L = 1. These yield to conclude that for z > a, the force 
difference is independent of z, such that 



AT 






V7r+^M 



x/3 



hu4 

V2 



= 0.202 



fojj„ - hio\ 



Hlu^ + hjj^ 



(25) 



proaches the plane. In the case of Al and Ag, the force 
difference at z/a = 0.001 is less than 0.01%, while at 
z/a = 0.1 the difference is about 1.2%. In Eq. Ij25(l we 
observe that as the plasma frequency of the materials 
becomes more dissimilar the difference of the force be- 
comes larger. For example, for potassium (hwp = 3.8 eV) 
and aluminum we would find that AT ~ 12.4% when 
2 > 2a, while for gold {ftw p = 8.55 eV) and cop- 
per (hujp = 10.8 eV), AT ~ 2.4% when z > 2a. In 
the other hand, the differences become very small when 
Z < 0.1 independently of the materials because of the 
high-multipolar surface plasmons interactions. 
V. SUMMARY 

In this work, we show that the Casimir force between 
a sphere and a plane can be calculated through the en- 
ergy obtained from their interacting surface plasmons. 
Our result is in agreement with the work of van Kam- 
pen and collaborators^, and GerlacbZ that showed that 
the Casimir force between parallel planes can be obtained 
from the zero-point energy of the interacting surface plas- 
mons of the planes. We show that the Casimir force of 
a sphere made of a given material A and a plane made 
of a given material B, is different from the case when the 
sphere is made of B, and the plane is made of A. We ob- 
tain a formula to estimate the force difference for Drude- 
like materials, when they are separated at least two radii 
of the sphere. We find that the difference depends on (i) 
the plasma frequency of the materials, and (ii) the dis- 
tance of separation between sphere and plane. We also 
found that as the sphere approaches the plane, the force 
difference becomes smaller. We conclude that the geom- 
etry of the system is important at distances larger than 
the radii of the sphere, where the force is dominated by 
the surface plasmons of the bodies. The energy of these 
surface plasmons is equal to the plasma frequency of the 
materials times a factor given by the geometry. 
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For the case of A=A1, and B=Ag, we find from Eq. I|25|) 
that AT ~ 4.9%, which is the calculated value shown 
in Fig. QJ) when z > 2a. When z < a the difference 
decreases as z/a also does, because more and more mul- 
tipolar interactions become important as the sphere ap- 
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APPENDIX 



The tensor that couples the interaction between the sphere and the plane i: 



,2:5 



A' 



I m 
Im 



Y^ m {0,v) 



i+i 



A+l'+i 



(47t) 3 (/ + 1 +m- m)\(l + 1 - m + m')\ 



{21 + l)(2l' + l)(2l + 21' + 1)(Z + m)\{l - m)\{l' + m')\(l' - m')\ 



1/2 



(A.26) 



where Yj m (0, <p) are the spherical harmonics. In spherical coordinates r= (2(z + a), 0, 0), then, the spherical harmonics 
reduce tc 31 



Y-r^) = 



2(/ + Q + l 
An 



1/2 



m—m ,0' 



(A.27) 



This means that m = m and multipolar moments with different azimuthal charge distribution are not able to couple 
or interact among them. This yields to 



Aim' 47r 

lm [2(z + a)] l + l '+ 1 



which is a symmetric matrix. 
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